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Abstract. Motivated by a proposal to create an optical helix-shaped waveguides 
for cold atoms and molecules, we discuss local perturbations which can create 
bound states in such a setting. This is known about a local slowdown of the twist; 
we show that a similar effect can result from a local tube protrusion or a change 
of the helix radius in correlation with its pitch angle. 



1 Introduction 

Geometrically induced bound states in tubular regions appear in various 
contexts - see, e.g., [DE951 fCDFK05j and references therein. In particular, 
effects of waveguide twisting were discussed recently in p^KK05. IEK05] . An 
independent strong motivation to look into such a problem came from a study 
[Bh07] which showed how helical waveguides for cold atoms and molecules 
can be created by means of two counter-propagating beams. Such waveguides 
do not have hard walls, of course, but it is known that spectrum-generating 
geometric mechanisms are robust - see, e.g., [EI01J - so it is useful to study 
such systems in a model which allows a mathematical treatment. 

It was observed in |EK05] that slowing down the waveguide twist acts 
as an effective attractive potential. Since the mentioned optical waveguides 
allow for various modifications it is useful to look for other perturbations 
which could have a similar effect. We will discuss two of them. One is a 
local protrusion which is well known to create bound states in straight tubes 
[BGRS97] . however, the proof needs to modified to fit to helical waveguides. 

Another possible perturbation is a variation of the helix radius. Here it is 
more suitable to employ the formalism introduced originally in [DE95] . An 
exhaustive analysis would be complicated, however, and we restrict here to 
the situation of a thin waveguide when the cross section size is much smaller 
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than the helix radius, paying attention to particular cases, those of a circular 
channel and of a flat "ribbon-shaped" one. In both cases we observe that 
a local reduction of the helix radius can induce a weakly bound state if the 
"pitch angle" is large enough, and in contrast, at small pitch angles localized 
states may appear for a local radius enlargement. At the same time, the 
crossover between the two regimes depends on the cross section. 

2 Helical channel Hamiltonians 

We will employ the notation of |EK05j which we extend to the case of varying 
cross-section and helix radius. Let u be an open connected set in R 2 we fix 
t° = (t°, t°) G to and define 

u(a) = {l a (t) := (a(t 2 - t° 2 ) + t° 2 , a(t 3 - t°) +4): t G u} ; (2.1) 

changing a means a radial scaling of the cross section u> w.r.t. the point t°. 
Let further 6(s), a(s) be differentiate functions from R to R. We define the 
mapping £ from 1 x w to I 3 by £ := C(s, t), where tilde! marks the scaling 
d2UD for a :=a(s), i.e. t= (t 2 , t 3 ) = (a(s)(t 2 - t° 2 ) + t° 2 , a(s)(t 3 - t° 3 ) + tf) , 
and 

C(s, t) := (s, t 2 cos6(s) +t 3 sin5(s), t 3 cos^(s) - i 3 sin 0(a)). (2.2) 

We denote the £-image of R x u by Q. It is tube in R 3 which is purely helical 
if the functions 9(s), a(s) are constant. We will be concerned with a quantum 
particle concerned within this tube assumimg that it has hard walls. Using 
appropriate units to get rid of unnecessary constants, we can then identify 
the particle Hamiltonian with the Dirichlet Laplacian H on L 2 (fl), i.e. the 
self-adjoint operator associated with the closed quadratic form 

Q[^] : = ! |VVfdtds, iieHliO). (2.3) 
Jn 

The family of operator^ corresponding to given 9(s), a(s) will be written as 
H(a). In particular, H (a) will denote the operator corresponding to 9(s) = 

would be appropriate, of course, to use t indexed by the scaling parameter and to 
mark its dependence on s. Such a notation, however would be rather cumbersome and we 
believe that no confusion will arise. The tilde will be later used for functions on the scaled 
cross section uj(a) as well. 

2 The function 9 will be kept fixed so we do not write it explicitly in the operator symbol. 
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fios, a(s), and H will refer to the purely helical tube, 9(s) = j3 s, a(s) = 
a . A simple substitution of variables shows that the quadratic form Qo(a) 
associated with H (a) is unitarily equivalent to 

Qo(a)M=/ / (|V^| 2 +|^ + ^:i 2 )dtd S , (2.4) 

JR iu(a(s)) 

denoted for simplicity by the same symbol, where we have introduced the 
notation f' T : = (t2<9* 3 — £2^)/ m accordance with [EK05J. 

Furthermore, the quadratic form (12.31) can transformed to other unitarily 
equivalent expressions supported by a straight tube (cylinder) is such a way 
that the geometric information is contained in the coefficients. Following 
[EKK05] we arrive at the quadratic form q' |EKK05] acting on L 2 (Q , IG] 1 ^ 2 ) 
as 

q'M:= / WfiG^id^lGl^dsdt, (2.5) 

where we number the variables in such a way that (<9i, <9 2 , <9 3 ) = (d s , d t2 , d tz ). 
By a straightforward computation we get for the metric tensor the expression 



/ \ h 2 h 3 \ 

a(s) a(s) 

Qij = 

V 
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(2.6) 



h 2 = (t 2 -t° 2 )a(s)+t 3 6(s), (2.7) 
h 3 = (t 3 -t° 3 )a(s)-t 2 0(s). (2.8) 

We can also pass to a quadratic form on the Hilbert space L 2 (R x oj) without 
the additional weight IG 1 ^ 2 ]; this is achieved by putting |EKK05| 

q[V] := q'[\G\-^} 

= {d?l>, G^d^) + (diF)G ij (djF)ijj) — 2Re(diip, G ij (d j F)iP)(2.9) 

where F := log(|(j| 1//4 ). If we finally plug in the metric tensor (12.61) we get 

qM = J — ^(K*)^ - h ^ ~ h ^ - «( s )^! 2 + iV^Hdsdt. 

(2.10) 
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After these preliminaries we shall formulate a theorem which will be 
proved in Section @] below. We will consider helical tubes with a constant 
"pitch angle", $(s) = (3qs for some (3q > 0. Furthermore, we assume that 
outside a compact region Q is purely helical; without loss of generality we 
may suppose that a(s) = 1 there so the unperturbed tube cross section is to. 

Theorem 1 Suppose that a(s) — l is a nonzero continuous function, which is 
nonnegative and compactly supported. Let further u, t° be such that uo{a) C 
u>(a') holds for a < a' , then the operator Ho (a) has at least one eigenvalue 
below the threshold of the essential spectrum. 

Remark We exclude here the case of a straight tube, /3q = 0. Bound states 
induced by a local protrusion exists in such a situation also, however, one 
can use a more straightforward way to prove the claim - cf . |BGRS97] . 



3 Spectrum of H 
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As in |EK05] our strategy is to regard H (a) as a perturbation of H = 
H (l). The spectrum of the latter operator is purely absolutely continuous 
and covers the half-line [-E'(l), oo), where E(a) is the lowest eigenvalue of 
the operator 

h(a) = -A^ {Q) - (3 2 (t 2 d ta - t 3 d t2 ) 2 (3.1) 

acting in L 2 (uj(a)) (pay attention to the fact that this h(l) corresponds to 
h(0) of |EK05j ). Moreover, the ground state is non-degenerate and the cor- 
responding ground-state eigenfunction f a {t) is strictly positive in u(a). 

A simple substitution of variables given by l a (12.11) allows us to pass to 
the unitarily equivalent operator h(a) acting on L 2 (uj), thus without a tilde. 
It acts as follows, 



(V, h(t*)V) = h(l)^) - 4 / 1*2^^ - hd t2 ^ 



_ |( a (t 2 _ t ° 2 ) + t° 2 )d t3 ^ - (a(t 3 - t° 3 ) + tl)d t2 ^\ 2 \ dt. (3.2) 

Next we ask about the dependence of h(a) on a. By f a (t) we denote the 
normalized ground-state eigenvector of h(a), i.e. f a (l a (t)) = f a {t). 
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Lemma 2 / a (£), E(a) are real- analytic functions of a in (0,oo). In particu- 
lar, there are < and such that E (a) = E + (a- + o(a- 1) 
and f a (t) = /(£) + (a - l)/«(t) + o(a - 1). 

Proof: We rewrite and estimate the integral on the right-hand side of ( 13. 2p . 



for some C\ > 0, since the last integral can be estimated by (diama; || V t ?/>||) 2 . 
It follows that the operators h(a) form an analytic family type of type (B) 
- cf. |Ka[ Sec. VII. 4], note that it is a particular case of exercise 4.23. 
there. Hence the analyticity of f a , E(a) follows from finite multiplicity of 
the ground state; recall that in fact it is non-degenerate. 

It remains to prove that E^ < which we will do using the minimax 
principle. Suppose that a > 1 and consider the test function obtained as a 
shifted ground state, ip(t) := f{t 2 + (a — l)a2, £3 + (a — 1)^3), where / means 
here the lowest eigenfunction of h(l) extended to u(a) by zero. For small 
enough shifts a 2 , (13 we have {£2 + (a — l)a 2 , £3 + (a — 1)03} C u(a) for all 
£ = (£1, £2) £ k>, and consequently 



If the term linear in (a — 1) does not vanish identically we are done. Suppose, 
on the contrary, that 




+ \a{t 2 d t3 iP - hd t2 i>) + (1 - a){t% 3 i; - t° 3 d t M 2 \ d£ 
<%-2a 2 | f \t 2 d t3 iP-t 3 d t2 iP\ 2 dt 




< Ci(V>, h(l)V) 



(V, h(a)V) = (/, h(l)/) 
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holds, then we employ another test function, for instance one obtained by 
scaling, $(t) := / + 1°, ^ + tjj) . Using JOD we get 

^,h(a)^) = E-f3 2 [ \t 2 d t .J-hd t J\ 2 dt-(l-\) [ \VJ\ 2 dt 
+ Pl I \t 2 d t J-t 3 d t jfdt + (a-l) 2 ^ I \t%J-t%J\ 2 dt 

J UJ ^ J UJ 

= E-2(a-l) [ |VJ| 2 dt + o(a-l), 
and since / \V t f\ 2 dt > we get the sought result. ■ 



4 Helical channel with a protrusion 

Now we are in position to prove Theorem [TJ By assumption the function 
a(s) — 1 is compactly supported, hence H (a) — H is relatively compact and 
Cess (Ho (a)) = Cess (Ho). Since we know the essential spectrum threshold, we 
can find eigenvalues below it using a variational estimate. 

Since a(s) — 1 is supposed to be nonzero and non-negative, one can find an 
interval (—Sq, Sq) within the support of this function on which the inequality 
1 + e\s — sq\ < a(s) holds for s small enough. It follows from the domain 
monotonicity of Dirichlet Laplacian that it is sufficient to establish existence 
of a bound state for a(s) := 1 + e\s — s |x(— s , s ). As usual in such cases 
we start constructing a trial function from then threshold-resonance of H . 
Given 5 > 0, e > we put ^ iE (s, t) = f E (t)4>s{s), where f e is the ground- 
state eigenfunction of h(l + e\s — so|x( — so 5 so)) and 

( e 5(s+So) if s < -s , 

(f> s (s) = I 1 if -s < s < s , 

[ e - s( - s ~ s ^ if s > s - 

We plug this expression into 1 12. 10p and by a straightforward computation we 
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get 



q[^s,e}-E\\^ 5 ^\\ 2 = 6+ / (E(l+e\ 8 -8 \)-E(l))ds 



2 Re 



-S 

so 1 A 







a(s 



{hdtje-kdtsfe 



x d s f £ 



a[s 



- 4)d t J £ + (t 8 - tg)flU - /.)) ^ dt 



+ 



so 



d s f £ - ^4 f & " ^R/s + (*3 - " /« 



dsdt. (4.1) 



Now we inspect the behavior of the expression for small e using Lemma [2j 
Using a(s) ~ esgn(s) and d s f £ ~ £Sgn(s)/W we obtain 



-left 



/SO 
E {1) \s~s \ds 
■so 



S() 



S JU 



X 



(/(i) _ {t2 _ _ (ts _ t o )dtJ _ ^ dt ds + o(£) _ (4 _ 2) 



The integral term vanishes due to the parity of the sign function; putting 
then 5 = e 2 we arrive at 

q[^ s J-E\\^ 6 , £ \\ 2 = sE^s 2 + o(e). 

omce is negative, it is sufficient to choose e small enough to conclude 
the proof. 



5 Thin helical tubes 

Next we look what happens if it is the radius rather than the cross section of 
the helical tube which is locally changed; for simplicity we restrict ourselves 
to small perturbations of a thin tube. To this aim it is more convenient to 
use the approach due to [GJ92] and [DE95] where the cross section is taken 
perpendicular to the tube axis. Let us stress that while helical tubes do not 
fall into the class of the asymptotically straight ones for which the existence 
of geometrically induced discrete spectrum was established in [D E95j . the 
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perturbation theory w.r.t. the tube radius developed there and in [UB96J 
remains nevertheless valid and we can use it here. 

The generating curve of our thin channel will thus helix a varying radius, 

f (t) = (t, R(t) cos 0(t), R(t) sin 9(t)) , tel, 

whose image in M. 3 will be denoted [r]. With F we conventionally associate 
its Frenet triad frame (t, n, b) consisting of its tangent, normal, and binor- 
mal vectors. Furthermore, k, t will denote the curvature and torsion of T, 
respectively. We suppose that 

9(t) = Pot 

and the radius is slightly changing according to 

R(t) = R + e5(t) , (5.1) 

where 5(t) is a fixed (nonzero) C 2 smooth function of compact support and 
e is a small parameter by which we mean that e||5||oo ^ Ro- 

For further reference, let us first inspect the unperturbed helix, e = ; 
the correspondent quantities will be indicated by the zero subscript. It is 
straightforward to check that for T (t) = (t, Rq cos j3 t, sin/5 t) the Frenet 
triad is 

1 RoPo sin p t R f3 cos f3 t 

cos fat, — sin Pot), 
Ropo sin p Q t cos Pot 




The normal vector no(t) is perpendicular to the helix axis while the tangent 
to(t) and binormal ho(t) contains with it nontrivial angles independent of t; 
the curvature and torsion are also constant and equal to 

RqPo Po 

K °-\+Ripr r °~i+m' 

To use the above mentioned results |DE95] and |CB96] we have replace t 
in the parametrization of [r] by the arc length of the curve, T(s) := T(t(s)), 
where t(s) is determined by the implicit equation 

rt{a) . 

|f(r)|dr. (5.2) 
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The parametrization change makes, of course, little difference for the unper- 
turbed helix where the two are mutually proportional, t(s) = s{\+I$q0q)~ x '' 2 . 

We will consider two models of thin helix quantum waveguides corre- 
sponding to different cross-sections. The latter will a family u(s) of bounded 
connected neighborhoods of zero, typically obtained by rotations of a fixed 
uj smooth w.r.t. s ; the tube in question Q C M 3 is then defined as the image 
of the map s i— > T(s) + x 2 n + x 3 b , (x 2 , x 3 ) G u(s) , as s runs through K. 
The Hamiltonian is again the Dirichlet Laplacian on Q denoted by H and 
the parametric description of Q makes it possible to replace it by a unitary 
equivalent operator on the "straightened" waveguide. The two models we 
will be interested in are the following: 

(i) The cross-section is two-dimensional, circular with T at its centre, and 
perpendicular to the helix. In that case we can make uo fixed in the 
so-called Tang coordinate system which rotates around t w.r.t. the 
Frenet triad with the angular velocity r. In that case one achieves a 
full decoupling of the longitudinal and transverse coordinates in the 
"straightening" transformation, see (DE95J for details. The perturba- 
tion theory with respect to the circle radius developed there shows, 
in particular, that the bottom of the spectrum for a thin tube is de- 
termined - after subtracting the continuum threshold energy - by the 
one-dimensional Hamiltonian T := — 4r + V^ TC (s) with the effective 
potential 

V^s) := -\k{s)\ (5.3) 

(ii) The optical waveguides which we use as a motivation [Bh07] are, how- 
ever, far of a circular shape having a very elongated cross section the 
sizes of which in two principal directions may differ by as much as 
two orders of magnitude. In such a case it is appropriate to use an 
idealized description due to [CB96J in which the cross section is a one- 
dimensional segment and Q has thus form of a winding ribbon; in accor- 
dance with [Bh07j we suppose that the segment u(s) is perpendicular 
to the helix axis. To achieve that, the function a describing the ro- 
tation of uo must be such that (ncosa — bsina)(s) is perpendicular 
the axis direction for any s, i.e. that the first component of this vec- 
tor vanishes. By the analysis of [CB96J the weak-coupling problem is 
again described in the leading order by a one one-dimensional operator 
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in which the effective potential (15.31) is replaced by 

V:t° n (s) := -\*\s) cos(a( S )) 2 + ~ (r(s) - a(s)f . (5.4) 

Conditions for the existence of weakly bound states can be thus deduced 
from well-known properties of effective operator — ^ + T4g which has isolated 
eigenvalues i^] the potential is attractive in the mean, J R (V e s (s) — E ) ds < 0, 
where E := lim| s |_ >00 V(s)); in our case the limit obviously exists since the 
curvature and torsion are constant outside a compact set. 

Let us look now what the above condition gives for the described geome- 
tries. Since the radius perturbation is weak by assumption, for the arc-length 
of perturbed helix we get from (15. 2p the relations 



1 + Rift t(s) + O(e) , t(s) = + 0(e) =: t (s) + 0(e) ; 

V 1 + H of j o 

since 5 G C 2 by assumption, also the first two derivatives of t(s) and to(s) 
coincide up to O(e). Then we can compute the geometric quantities which 
enter the above expressions for effective potentials; after a straightforward if 
tedious computation we get 



, c x , (/3 2 -«)^o( g ))-K/j 2 + l)^ ( g )) ^ 

( ' s) " + (i + RlPiy £ + U{£ } 

o iZg^(t (a)) + (R 2 (3 2 + l)S(t (s)) ^ 2i 
T(S) " T ° " 2 RMl + RlPl? £ + ° {£l 



tan a{s) = 6{t ° {s)) e + 0(e 2 ) 



Comparing this results with the effective potential for the ribbon ( 15. 41) we 
see that the terms linear in 8(to(s)) do not contribute to the integral in 
the condition J(V e s(s) — E ) ds — if 8 is smooth as assumed, and thus 
the ribbon twisting described by the function a plays no role in the leading 
order. Computing the effective potentials explicitly with the help of the 
above formulae we arrive at the following conclusions: 



3 We leave out the critical case, J (V c g(s) — Eq) ds = 0, since using the effective operator 
we deal the leading order only. 
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(i) for the circular tube we have 



eff 



4(1 + ^2)2 



Ro^jRlPj - l)6(t(s)) + RoP 2 (R 2 P 2 + l)5(t(s)) , 
+ 2(1 + Rift)* e + U(s). (5.5) 

Hence we can distinguish two cases: for a "steep" helix, RqPo > 1, a 
weakly bound state occurs if L S(t) dt < , i.e. in the situation where 
the the helix radius is locally reduced. On the other hand, for a small 
pitch angle, Ro/3 < 1 the bound state occurs if j R S(t)dt > 0, i.e. if 
the radius is locally enhanced. 

(ii) in the ribbon case we have 

T/ribbon^ _ R-ofto — 2/?q 



" ir V; " 4(1 + J R^ )2 

- 5)*(t(*)) + (1 + ^/9o 2 )(^/^o 2 - 2)^(t(a)) 



+ 



2 J R (l + ^/?o 2 ) 3 

+ 0(£ 2 ), (5.6) 



and again we have two cases differing from the previous situation just 
by the critical value of the pitch angle. For RoPo > V5 the bound state 
occurs under the local "squeezing", while for RqPo < Vo we have to 
"inflate" the helix locally to achieve binding. 
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